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Abstract:  Sensitivity  analysis  in  data  envelopment  analysis  for  the  Chames-Cooper- 

Rhodes  ratio  model  is  studied  for  the  case  of  the  simultaneous  proportionate  increase  of 
all  inputs  and  proportionate  decrease  of  all  outputs  of  an  efficient  decision  making  unit 
for  which  efficiency  is  preserved.  Sufficient  conditions  which  preserve  efficiency  are  found 
and  a  numerical  example  illustrating  the  results  is  provided. 
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1.  Introduction 

Sensitivity  analysis  in  Data  Envelopment  Analysis  (DEA)  for  the  Charnes- 
Cooper- Rhodes  (CCR)  ratio  model  was  introduced  by  Charnes  et  al.  [5]  for  the  case 
of  the  change  of  single  output.  Sufficient  conditions  for  an  efficient  Decision  Making 
Unit  (DMU)  to  continue  to  be  efficient  after  the  change  of  single  output  were  found 

*This  research  was  partly  supported  by  NSF  Grant  SES8408134and  ONR  Contract  N00014-82- 
K-0295  with  the  Center  for  Cybernetic  Studies,  The  University  of  Texas  at  Austin.  Reproduction 
in  whole  or  in  part  is  permitted  for  any  purpose  of  the  United  States  Government.  This  research 
was  also  partially  supported  by  the  Research  Council  of  Croatia,  Yugoslavia. 

The  paper  is  prepared  for  presentation  at  1 1th  European  Congress  on  Operational  Research 
EURO  -  XI,  July  16  -  19,  1991,  Aachen,  Germany. 
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first  [5].  The  generalizations  of  that  result  for  the  case  of  the  simultaneous  change 
of  all  outputs,  the  case  of  the  simultaneous  single  output  and  single  input  changes, 
the  case  of  the  simultaneous  change  of  all  inputs  and  the  case  of  the  simultaneous 
change  of  all  inputs  and  outputs  for  the  CCR  ratio  model  were  given  by  Charnes 
and  Neralic  [6],  [7],  [8].  Similar  results  for  the  additive  model  were  found  by  Charnes 
and  Neralic  [9].  Sufficient  conditions  for  an  efficient  DMU  to  preserve  efficiency  after 
the  proportionate  change  of  inputs  (or  outputs)  were  giwen  by  Charnes  and  Neralic 
[10].  The  results  in  the  case  of  the  proportionate  change  of  inputs  (or  outputs) 
can  be  used  for  ranking  among  efficient  DMUs  as  it  was  suggested  by  Banker  and 
Gifford  [2]. 

The  aim  of  that  paper  is  to  study  the  case  of  the  simultaneous  proportionate 
change  (increase)  of  inputs  and  proportionate  change  (decrease)  of  outputs  of  an 
efficient  DMU  preserving  efficiency.  Using  the  results  of  Charnes  and  Neralic  [8] 
in  sensitivity  analysis  in  DEA  for  the  CCR  ratio  model  sufficient  conditions  for  an 
efficient  DMU  to  preserve  efficiency  after  the  simultaneous  proportionate  change  of 
inputs  and  outputs  are' given. 

The  paper  is  organized  as  follows.  The  results  in  sensitivity  analysis  in  DEA 
which  will  be  used  later  are  contained  in  Section  2.  The  main  result  of  the  paper  is 
given  in  the  Theorem  2  in  Section  3.  Section  4  gives  an  illustrative  example.  The 
last  Section  contains  some  conclusions. 


2.  Preliminaries 

Let  us  suppose  that  there  are  n  Decision  Making  Units  (DMUs)  with  m  inputs 
and  s  outputs.  Let  xtJ  be  the  observed  amount  of  zth  type  of  input  of  the  jth  DMU 
(  Xij  >  0,  i  =  l,2,...,m,j  =  1,2,  ...,n)  and  let  yrj  be  the  observed  amount  of 
output  of  the  rth  type  for  the  jth  DMU  (yrj  >  0,  r  =  1,2, . . .  ,s,j  =  1,2, .. .  ,n). 
Let  Yj,  Xj  be  the  observed  vectors  of  outputs  and  inputs  of  the  DMU,,  respectively, 
j  =  1,2, ...  ,n.  Let  e  be  the  column  vector  of  ones  and  let  T  as  a  superscript  denote 
the  transpose.  In  order  to  see  if  the  DMUjo  =  DMUo  is  efficient  according  to  the 
CCR  ratio  model  the  following  linear  programming  problem  should  be  solved: 

min  0At  +  •  •  •  +  0Ao  +  •  •  •  +  0An  —  eeTs+  —  eeTs~  +  0 

subject  to 

V'iAi  -f  •••  +  VoAo  +  •••  +  Yn An  —  s+  =  Vo 

— A  i  A  |  —  —  A0A0  —  ...  —  Xn\n  —  s  +  X0  0  —  0  (1) 

A),  -  •  - ,  A„,  s  ,  s  ^  0, 

with  Y0  =  Y]0 ,  X0  =  XJg.\0  =  A ]0  and  0  unconstrained.  The  symbol  t  represents 
the  infinitesimal  we  use  to  generate  the  non- Archimedean  ordered  extension  field  we 


2 


shall  use.  In  this  extension  field  t  is  less  than  every  positive  number  in  our  base  field, 
but  greater  than  zero.  DMU0  is  DEA  efficient  if  and  only  if  for  the  optimal  solution 
(A*,s+’,s~*,0*)  of  the  linear  programming  problem  (1)  both  of  the  following  are 
satisfied  (for  details  see  [3]): 

min  0  =  O'  =  1 

.s+*  =  s~m  =  0,  in  all  alternative  optima.  ' 

We  are  interested  in  variations  of  all  inputs  and  all  outputs  of  an  efficient 
DMUo  preserving  efficiency  A  decrease  of  any  input  cannot  worsen  an  already 
achieved  efficiency  rating.  Downward  variations  of  inputs  are  not  possible  in  the 
efficiency  rating  for  an  efficient  DMUo  *  Hence  we  can  restrict  attention  to  upward 
variations  of  inputs  of  an  efficient  DMUo  which  can  be  written  as 

F«o  =  x,o  +  ft,  ft  >  0,  i  =  1,2, ...  ,m.  (3) 

Similarly,  an  increase  of  any  output  cannot  worsen  an  already  achieved  efficiency 
rating.  Upward  variations  of  outputs  are  not  possible  in  the  efficiency  rating  for  an 
efficient  DMUo-  Hence  we  can  restrict  attention  to  downward  variations  of  outputs 
which  can  be  written  as 


Vro  -  !/ro  -  ar  >  0,  aT  >  0,  r  =  l,2,...,s.  (4) 

For  an  efficient  DMU0  because  of  (2)  vectors  [  Y0  -  X0  }T  and  [0  X0  ]T  must 
occur  in  some  optimal  basis,  which  means  that  there  is  a  basic  optimal  solution 
to  (1)  with  Aq  =  1  and  0‘  =  1.  Changes  (3)  and  (4)  are  then  accompanied  by 
alterations  in  the  inverse  B~}  of  the  optimal  basis  matrix 


Yb  -/£  0  0  ' 

~XB  0  -Is  XQ  ’ 


(5) 


which  corresponds  to  the  optimal  solution  (A*,  s+*,  s~",  0")  of  (1)  with  Aj  =  1  and 
0*  =  1 .  Let 

»  i,  j  -  l,2,...,s  +  m, 

be  the  inverse  of  the  optimal  basis  B  in  (5).  Let  P},  j  =  1,2, ..., n  +  s  +  m  + 1  be  the 
columns  of  the  matrix  and  let  P0  be  the  right  hand  side  of  the  linear  programming 
problem  (l).We  will  use  the  following  notations: 


=  B  lPj,  j  =  0, 1,. ..  ,n  +  5  +  m  +  1, 
u.'r  =  ctbB~\ 

=  ctbB-'P3 

=  u1  Pj,  j  =  0, 1,. . .  ,n  +  s  +  m  +  1. 
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The  simultaneous  change  of  outputs  (3)  and  inputs  (4)  leads  to  the  following  change 
of  the  optimal  basis  matrix  B 

B  =  B  +  AB  (6) 


A£  = 


0 

...  o 

-<*1 

0  ... 

0 

0 

...  0 

~a2 

0 

0 

0 

...  o 

~0ts 

0  ••• 

0 

0 

...  o 

-01 

0  ••• 

01 

0 

...  o 

-02 

0  ••• 

02 

0 

...  o 

-0m 

0  ... 

0m 

and  the  following  change  of  the  right  hand  side  vector 

Bo  —  Bo  +  [-£*1  —  Qf2  . . .  —  as  0  . . .  0  ]T  , 

where  indexes  k  and  s  +  m  correspond  to  the  optimal  basic  variables  AJJ 
O'  =  1  respectively.  Using  matrices 

<*i 

02  O2 


(8) 
1  and 


(5+m)x2 


o,  a, 

0i  0 
02  0 


,0m  0 


k  i+m 

VT  _  0  0  —  1  0  •••  0  1  /in\ 

2x  (s+tn)  [o  •  •  0  0  0  •••  0  -1  '  ’ 

we  can  write  the  perturbation  matrix  (7)  as  A B  =  UVT .  Let  us  use  the  abbreviation 

M  =  I  +  VTB~lU , 

where  matrix  M  is  nonsingular  with 

5  m 

del  M  =  1  -  +  K+m,s+t)0t  + 

t=i  «=i 

+  (±K:m,ot)(±KUA)  -  (iK^t)(±b:lm^t0t),  (11) 
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and 


D  =  UM~1VT. 


(12) 


Theorem  1.  Conditions 

uT DY j  >  Zj  -  Cj,  j  an  index  of  nonbasic  variables ,  (13) 

are  sufficient  for  DMU0  to  be  efficient  after  the  simultaneous  changes  of  inputs  (3) 
and  of  outputs  (4).  If  detM  >  0,  conditions  (13)  can  be  written  in  the  following 
way 

7*  IV,  +  7a+mrj+miJ  >  (zj  -  Cj)det  M,  (14) 

with 

ms  am 

3k  =  “(I  +  +  I"1  +  XIX+m.tat)(XIUW?t)>  (15) 

(=1  (=1  «=1  t=l 

and 

ms  am 

7  ,+m  =  (Y^bk!s+A)(J2u^t)  +  (1  ~^2bf}at)(J2^+t/3t).  (16) 

(=1  <=1  t=l  t= 1 

For  the  proof  and  details  see  [8]. 


3.  Simultaneous  proportionate  change  of  inputs 

and  outputs 

Let  us  consider  the  simultaneous  proportionate  change  (increase)  of  all  inputs 

•r,o  =  0x,o,  0>l,  i  =  l,2,...,m,  (17) 

and  the  proportionate  change  (decrease)  of  all  outputs 

g^ o  cxyrQ,  0  Q  ^  1,  r  —  1,2, ...  ,s,  (16) 

of  an  efficient  DMUo  preserving  efficiency.  We  are  interested  in  sufficient  conditions 
for  DMU0  to  preserve  efficiency  after  the  simultaneous  changes  (17)  and  (18). 

Theorem  2.  .  Let  us  suppose  that  DMU0  is  efficient  and  let 

det  M  =  1—  <zj(l  —  or)  +  (—  61+62 )(/?  —  1)  +  ( a2^i—  ai^2)(l—a)(^~l)^>0,  (19) 

with 

s  s  mm 

=  Yl^ktUlO-  «2  =  bl  =  XI  Vl+fX‘°’  ^2  =  XI  ^3+m,s+tXtO-  (20) 

t=l  (=1  t= 1  t=l 
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a3  —  0,  ^3  —  (21) 

£=1  t=l 

dj  =  — <2.3  r  i-j  +  aiCj,  t}  =  —  b3(Tkj  —  ri+mj)  +  (— bi  +  &2  )c;>  (22) 

f:  —  (a2^3  —  “3^2  JTfcj  +  (03^1  —  ai^3)r3+mj  —  ( CI2&1  —  ai^2)Cj,  (23) 

J  =  1,2,.. .,n  +  5  +  m  +  1, 

with  Cj  —  Zj  —  Cj.  Then  the  conditions 

d}(  1  -  5)  +  Cj(0  -  1)  +  £(1  -  5)0  -  1)  >  Cj,  (24) 

j  an  index  of  nonbasic  variables , 


are  sufficient  for  DM Uq  to  preserve  efficiency  after  the  simultaneous  proportionate 
changes  of  inputs  (17)  and  of  outputs  (18). 

Proof: First  of  all  let  us  show  that  the  proportionate  changes  (17)  and  (18)  are 
the  special  cases  of  the  changes  (3)  and  (4)  respectively.  Using  the  substitutions 

fl=l+0,  0>O,  (25) 

and 

J,  =  d.r,o,  0i>  0,  i  —  1,2, ...  ,m,  (26) 

we  can  write  (17)  as 


■?,0  —  £|0  +  P^iO 

=  1,0  +  /?,,  Pi  >  0,  1  =  1,2, 


(27) 


It  means  that  the  proportionate  change  of  inputs  (17)  is  the  special  case  of  the 
change  of  inputs  (3)  with  p,,i  =  1,2,  ...,m  in  (26)  and  0  in  (25).  Similarly,  if  we 
put 

q  =  1  —  a,  0  <  a  <  1,  (28) 


and 


Qr  =  «J/ro,  Q  >  0,  t  =  1,2,  ...  ,5, 


(29) 


we  can  write  (18)  as 


VtO  =  2/rO  -  «2/v 0 

=  Uro  ~  Gr  >  0,  Qr  >  0,  r  —  1,2, ...  ,s.  (30) 

It  means  that  the  proportionate  change  of  outputs  (18)  is  the  special  case  of  the 
change  of  outputs  (4),  with  ar,r  =  1,2,. . .  ,s  in  (29)  and  a  in  (28). 
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Let  us  suppose  that  conditions  (24)  are  satisfied.  Then  using  (25),  (28),  (20)- 
(23)  and  (19)  it  is  easy  to  show  that  conditions  (24)  are  equivalent  to  conditions  (14) 
for  the  case  with  /?,,i  =  1,2, ...  ,m  in  (26)  and  ar,r  =  1,2,. . .  ,s  in  (29).  According 
to  Theorem  1  conditions  (14)  are  sufficient  for  DMUo  to  preserve  efficiency  after  the 
changes  (3)  and  (4).  Because  of  the  equivalency  between  conditions  (14)  and  (24) 
for  the  special  case  with  B, ,  i  =  1,2,  ...,m  in  (26)  and  Q;r,r  =  1,2,  in  (29), 
which  means  the  simultaneous  proportionate  changes  of  inputs  (17)  and  outputs 
(18),  it  follows  that  conditions  (24)  are  sufficient  for  DMUo  to  preserve  efficiency 
after  the  simultaneous  proportionate  changes  of  inputs  (17)  and  outputs  (18)  and 
completes  the  proof. 

Remark  1.  For  the  case  det  M  <  0  instead  of  det  M  >  0  in  (19),  the  inequality 
sign  >  in  conditions  (24)  should  be  changed  into  <. 

Remark  2.The  system  of  inequalities  (24)  together  with  conditions  (17),  (18) 
and  (19)  for  a  and  @  gives  the  area  AQ  in  the  plane  with  the  coordinate  system 
aO/?.For  each  point  (q,  /3)  in  the  area  A0  efficiency  of  DMUo  will  be  preserved  after 
the  simultaneous  proportionate  changes  of  inputs  (17)  and  outputs  (18). 

Remark  3.  We  can  use  the  area  A0  for  ranking  among  efficient  DMUs.  For 
example,  if  for  efficient  DMU|  and  DMU2  holds  A\  >  A2  it  can  be  said  that  ”DMUi 
is  relatively  more  efficient  than  DMUj”  because  DMUi  is  less  sensitive  to  the  si¬ 
multaneous  proportionate  change  of  inputs  and  outputs  preserving  efficiency  than 
DMU?.  The  ranking  among  efficient  DMUs  can  also  be  based  on  the  proportionate 
change  of  inputs  (or  outputs)  as  it  was  suggested  by  Banker  and  Gifford  [2]  and 
used  by  Charnes  and  Ncralic  [10]. 

4.  Illustrative  example 

We  will  consider  the  following  example  taken  from  [11]  with  five  DMUs,  one 
output,  two  inputs  and  data  in  Table  1. 

Table  1 

Data  for  the  example 


DMUj  1  2  3  4  5 

Output/Input 

^  2  4  2  3  2 

4  12  8  6  2 

.v  2j  6  8  2  6  8 

We  are  interested  in  the  efficency  of  DMU<,  with  Xq  =  [6  6]  and  Yo  =  [3].  In 
order  to  see  if  DMU4  is  efficient,  the  following  linear  programming  problem  should 
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be  solved: 


subject  to 


minOAi  4  OA2  +  OA3  4  OA4  4  OA5  —  £s+  —  £s4  —  £s2  4  6 


2Ai  4 

4  A2 

4 

2A3 

4 

3A4  4 

2A5 

— 

e  + 

S1 

= 

3 

— 4Aj  - 

12A2 

— 

8A3 

— 

6 A.,  - 

2A5 

-  -r 

46  0  = 

0 

— 6AX  - 

sa2 

- 

•2A3 

— 

6A4  - 

8A5 

- 

—  s2  4  60  ~ 

0 

A 1 , 

A2! 

A31 

A4, 

Asi  st , 

*r, 

H 

>  0. 

The  optimal  solution  of  problem  (31)  is 

Aq  = 

■K 

=  1 ,0-  = 

1,  a;  =  a;  = 

A: 

(31) 


0,  s|*  =  si’  =  s2’  —  0  with  min 0  =  O'  =  1,  which  means  that  DMXJ4  is  efficient. 
The  optimal  basic  variables  are  A3,  A4  and  8.  The  optimal  basis  matrix  is 


B  = 


2  3  0 

-8  -6  6 
-2  -6  6 


with  the  inverse 


/r1  = 

and  corresponding  optimum  tableau  in  Table  2. 


r  0  -I 

I  1 
3  9 

I  j_ 
3  18 


I  1 
6 
1 
9 
I 
9 


(32) 


Table  2 

Optimum  tableau 


r, 

r2 

1 3 

r5 

r6 

r7 

r8 

t9 

r0 

a3 

I 

3 

2 

3 

1 

0 

-1 

0 

1 

6 

1 

6 

0 

0 

8 

9 

8 

9 

0 

1 

1 

3 

1 

3 

1 

9 

1 

9 

0 

1 

0 

7 

9 

2 

9 

0 

0 

1 

_  3 

1 

5 

1 

18 

1 

9 

1 

1 

Zj  —  Cj 

2 

_  2 

9 

0 

0 

1 

-  3 

T  3  +f 

-Z^±£- 

~9  +  £ 

0 

1 

Let  us  consider  the  simultaneous  proportionate  change  (increase)  of  inputs 

•?io  =  60,  .r20  =  6)0,  0  >  1  (33) 

and  proportionate  change  (decrease)  of  output 

yio  =  3a,  0  <  a  <  1,  (34) 

of  DMU4  preserving  efficiency.  Using  (25)  -  (26)  in  (33)  we  get 

*io  =  640i,  0i=60,  0i  >0,  (35) 
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(36) 


i^O  —  6  +  02,  02  ~  60,  02  >  0. 

Similarly  using  (28)-(29)  in  (34)  we  get 

t/10  =  3  -  a,  >  0,  a,  =  3a,  0  <  a,  <  3.  (37) 

Using  (32)  we  have 

U >T  =  CT8B~ 1  =  [0  0  1]  B~l  =  [1/3  1/18  1/9].  (38) 

If  we  use  (32),  (38),  s  =  1,  m  =  2,  fc  =  2,  s  +  m  =  3  and  the  elements  of  Table  1 
it  is  easy  to  get 

a,  =  1,  «2  =  1,  &i  =  0,  b2  =  1,  a3  =  1,  63  =  1.  (39) 

Because  of  (39)  it  follows  from  (19),  (22)  and  (23) 

dct  M  =  1  -  (1  -0)  +  (0-  1)  -  (1  —  a)(0  —  1)  >  0,  (40) 

(i]  =  —  r2j  +  cj ,  ej  =  — (^2;  ~  r3j)  —  Cj ,  (41) 

and 

/,  =  -r3,+c„  _/  =  1,2, ...  ,9,  (42) 

respectively,  with  c}  =  —  c^.  Using  (41)  and  (42)  the  conditions  (24)  can  be 

written  in  the  following  way 

(-r2j+ci)(i -a)  +  (-r2j  +  r3j-cJ)(^-i)  +  (-r3j+cJ)(i-a)(^-i)  >  cJ?  (43) 

J  =  1,2, 5, 6, 7, 8. 

For  example,  if  j  =  1  using  elements  of  Table  2  we  have  from  (43) 

(-8/9  -  2/9)(  1  -  o )  +  ( —8/9  -  2/9  +  2/V)(0-  1 )  +  (2/9  -  2/9)(l  - &){0 -  1)  >  -2/9, 
or 

0  <  1.255. 

It  is  easy  to  see  that  the  solution  set  of  the  system  of  inequalities  (43)  together  with 
constraints  (33),  (34)  and  (40)  is  the  triangle  ABC  in  the  plane  with  the  coordinate 
system  5 00  in  Figure  l,  with  A(l,  1.25), 

Figure  1  about  here 


9 


B(0.8,  1)  and  C(  1,  1).  The  redundant  constraints  are  not  sketched  in  Figure  1.  For 
every  point  {a,0)  which  belongs  to  the  triangle  ABC  the  efficiency  of  DMU<  will 
be  preserved  after  the  simultaneous  proportionate  change  of  inputs  (33)  with  the 
coefficient  0  and  proportionate  change  of  output  (34)  with  the  coefficient  a. The 
point  C(l,  1)  means  that  there  are  no  changes  of  inputs  and  of  output,  the  point 
A(l,  1.25)  means  the  maximal  proportionate  increase  of  inputs  of  DMU4  for  25% 
preserving  its  efficiency  and  the  point  B(0.8,  1)  means  the  maximal  proportionate 
decrease  of  output  of  DMU^  for  20%  preserving  efficiency  of  DMU*.  These  results 
of  proportionate  change  of  inputs  (or  output)  in  that  example  are  the  same  as  in 
Charnes  and  Neralic  [10],  but  as  can  be  seen  in  Figure  1  these  changes  can  not  be 
done  simultaneously. 

The  area  X4  =  0.025  of  the  triangle  ABC  can  be  used  for  ranking  DMU* 
among  the  other  efficient  DMUs.  We  can  consider  DMU3  and  DMU5  which  are 
efficient  too.  It  is  easy  to  show  that  in  the  case  of  the  simultaneous  proportionate 
change  of  inputs  and  output  of  DMU3  preserving  efficiency  for  the  corresponding 
area  holds  /13  =  0. 19082  (for  -  =  0.00001).  In  the  same  way  it  easy  to  seee  that  for 
the  efficient  DMU5  holds  /Vs  =  0.49982  (for  e  =  0.00001).  According  to  the  Remark 
3  because  of  A3  =  /Vs  >  /l.t  it  means  that  ”DUM3  and  DMU5  are  relatively  more 
efficient  than 


5.  Conclusions 

The  simultaneous  proportionate  change  of  inputs  and  proportionate  change 
of  outputs  of  an  efficient  DMU0  preserving  efficiency  in  the  case  of  the  CCR  ratio 
model  in  DFA  is  studied  in  the  paper.  Using  the  results  of  Charnes  and  Neralic  [8] 
in  sensitivity  analysis  in  DFA  for  the  CCR  ratio  model  sufficient  conditions  for  an 
efficient  DMU0  to  preserve  efficiency  are  established  for  the  case  of  the  simultaneous 
proportionate  increase  of  inputs  and  proportionate  decrease  of  outputs.  Sufficiency 
conditions  give  lor  each  efficient  DMUo  the  area  which  can  be  used  for  ranking 
among  efficient  DMUs.  A  numerical  example  illustrating  the  results  is  provided. 

The  simultaneous  proportionate  change  of  inputs  with  the  coefficient  0  and 
proportionate  change  of  outputs  with  the  coefficient  a  which  is  studied  can  be 
generalized.  For  example,  the  cases  of  the  proportionate  change  of  inputs  with 
different  coefficients  0,.  1  =  1,2, ...  ,m  or/and  the  proportionate  change  of  outputs 
with  different  coefficients  or,r  =  1,2, ...  ,s  can  be  considered.  These  cases  seems  to 
be  interesting  also  for  the  BCC  model  [1]  and  the  additive  model  [4]  .  The  results 
for  these  cases  will  be  presented  elsewhere. 
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